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Abstract 

We use early-time observations of He I 10830 A to measure the extent of up- 
ward mixing of radioactive material in SN 1987A. This work develops and extends 
the work of Graham (1988), and places constraints on actual explosion models. 
The presence of the He I 10830 A (2s 3 S-2p 3 P) line at > 10 days post-explosion 
implies re-ionisation by 7-rays from upwardly-mixed radioactive material pro- 
duced during the explosion. Using the unmixed explosion model 10H (Woosley 
1988) as well as mixed versions of it, we estimated the 7-ray energy deposition 
by applying a purely absorptive radiative transfer calculation. The deposition 
energy was used to find the ionisation balance as a function of radius, and hence 
the 2s 3 S population density profile. This was then applied to a spectral synthesis 
model and the synthetic spectra were compared with the observations. Neither 
model 10H nor the mixed version, 10HMM, succeeded in reproducing the ob- 
served He I 10830 A line. The discrepancy with the data found for 10HMM is 
particularly significant, as this model has successfully reproduced the X-ray and 
7-ray observations and the UVOIR light curve. We find that a match to the 
He I line profile is achieved by reducing the extent of mixing in 10HMM. Our 
reduced-mixing models also reproduce the observed 7-ray line light curves and 
the iron-group velocities deduced from late-time infrared line profiles. We sug- 
gest that the He I line method provides a more sensitive measure of the extent 
of mixing in a type II supernova explosion. 
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1 Introduction 



It has long been recognised (Falk & Arnett 1973) that hydrodynamic instabilities should 
occur in type II supernova explosions. Chevalier (1976) studied the stability of shocks 
propagating through simple power-law density distributions and showed that for some 
cases the matter behind the shock is subject to instabilities. Bandiera (1984) further 
emphasised the role of such instabilities in generating chemical mixing. 

SN 1987A provided the first observational evidence for extensive mixing of the ejecta. 
Simple stratified models did not account for the extended, smooth, rounded plateau in 
the UVOIR light curve which was observed during the era covering 25 to 125 days post- 
explosion. To successfully model the light curve, it proved necessary to invoke upward 
mixing of radioactive nickel and cobalt to high velocities together with downward mixing 
of hydrogen deep into the core (Arnett 1988, Woosley 1988, Shigeyama et al. 1988, 
Shigeyama & Nomoto 1990, Nomoto et al. 1998). Upward mixing of radioactive 56 Ni 
and 56 Co to the outer regions of the envelope was also indicated by the early detection 
of X-rays (e.g. Dotani et al. 1987, Sunyaev et al. 1987) and 7-rays (e.g. Matz et al. 
1988, Mahoney et al. 1988, Sandie et al. 1988, Cook et al. 1988). Furthermore, the 
expansion velocities inferred from the line widths of infrared spectral lines of Fe II and 
Co II (Spyromilio, Meikle & Allen 1990, Haas et al. 1990, Tueller et al. 1990) indicated 
that a fraction of the iron-group elements had been dredged up to velocities as high as 
-3000 km/s. 

These observations stimulated the modelling of explosions both in two dimensions (Ar- 
nett et al. 1989, Hachisu et al. 1990, Fryxell et al. 1991) and three dimensions (Muller 
et al. 1991). This theoretical work showed that significant mixing occurs as the ex- 
plosion shock plows through the ejecta. However, the maximum iron-group velocities 
derived did not exceed 1300 km/s. This discrepancy with the observations was at- 
tributed to non-inclusion in the models of the radioactive decay energy deposited in 
the ejecta. Herant & Benz (1991) included the radioactive decay energy in their two- 
dimensional code, resulting in the formation of a giant "nickel bubble" during the first 
few weeks. In spite of this, their model still only attained iron-group velocities up 
to ~2000 km/s - 1000 km/s less than that indicated by the observed iron-group line 
profiles. 

Herant & Benz (1992) excluded the possibilities that the velocity discrepancy was due 
to three-dimensional effects, errors in the adopted progenitor structure, or numerical 
resolution. They argued that in order to explain the observed high velocities, it was 
necessary to invoke mixing during the early stages of the explosion when nickel is created 
by explosive nucleosynthesis (Herant et al. 1994). By pre-mixing nickel out to 1.5 M 
above the mass cut, they found that the subsequent Rayleigh- Taylor instabilities could 
produce iron-group velocities as high as those observed. Such pre-mixing is possible 
within the context of the delayed explosion mechanism (Herant, Benz, & Colgate (1992); 
Herant et al. 1994). In this scenario, the explosion is driven by neutrinos from the 
proto-neutron star which deposit energy in material just above the neutrinosphere, 
producing a hot, high-entropy bubble interior to low-entropy shocked matter. This 
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situation is prone to strong convective motions under the influence of the gravitational 
pull of the proto-neutron star, thus breaking the spherical symmetry. The residual 
inhomogeneities due to these early instabilities may be sufficient to seed subsequent 
shell-interface Rayleigh- Taylor instabilities. 

Another way to account for the observed high velocities is through the effect of a 
spherically symmetric shock on irregularities or "seeds" formed in the 16 burning 
shell before core-collapse. Recent hydrodynamic calculations (Bazan & Arnett 1998, 
Arnett 1994) have shown that the oxygen burning is intermittent, chaotic and strongly 
localised. The heterogeneous composition developed provides "seed" perturbations in 
density. These are believed to be sufficiently large to produce hydrodynamic instabilities 
which cause mixing in the presupernova envelope. This occurs in precisely the region 
where 56 Ni is explosively produced by oxygen burning behind the explosion shock. 
However, quantitative results are not yet available. 

To test the possible mixing mechanisms described above we need to be able to probe the 
conditions deep inside the supernova. In principle, much can be learned from the 7-ray 
and X-ray light curves. In addition, once the ejecta have become optically thin (t> 150 
days), high signal-to-noise infrared spectral line profiles of heavy elements can provide 
a powerful probe of their spatial distribution, and degree of mixing (Spyromilio, Meikle 
& Allen 1990, Herant & Woosley 1994). Unfortunately, apart from the exceptionally 
close SN 1987A, current instruments are insufficiently sensitive to allow such studies to 
be made for typical nearby type II supernovae. 

An interesting alternative is to probe the mixing of the nickel using infrared spectra at 
early times. This approach was first presented by Graham (1988), and was developed 
further by us (Fassia et al. 1998). The method is based on the detection and modelling 
of the He I 10830 A (2s 3 S-2p 3 P) line. This is a high excitation line. It is formed in a 
region where the population of the metastable 2s 3 S level is maintained by the balance of 
the recombination rate of He + with the rates of collisional de-excitation and forbidden 
radiative decay. However, because the He + recombination timescale is < 1 d (Graham 
1988) we expect that at late times (t > 10 d) all helium should be neutral and in its 
ground state. Thus, if He I 10830 A is present at times after 10 d there must exist 
a re-ionising mechanism. The identification of this line well after 10 d in SNe 1987A 
and 1995V led Graham (1988) and Fassia et al. (1998) to propose that the ionisation 
was maintained by the 7-rays emitted from the decay of 56 Ni/ 56 Co. The subsequent 
recombination produces the He I 10830 A line. Owing to the exponential sensitivity 
of the 7-ray flux to the column depth (Pinto & Woosley 1988) the observed He I line 
profile is sensitive to the amount of radioactive material moving at high velocities. 
Thus, the He I 10830 A line can probe the extent of upward mixing or "dredge-up" 
of iron-group elements in a type II supernova. This was demonstrated by Fassia et al. 
(1998) for the type Hp supernova SN 1995V. While an unmixed explosion model did 
not reproduce the observed He I line, by invoking upward mixing of about 10 -6 M of 
56 Ni to velocities above 4000 km/s (i.e. above the helium photosphere) we were able to 
match the observed spectra. 

In the present work we use the He I line method to carry out a detailed study of the 
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dredge-up of radioactive material in the ejecta of SN 1987A, and to place constraints on 
possible mixing models. The method is described in section 2. In section 3 we present 
results from the comparison of theoretical models with the data and estimate the 56 Co 
dredge-up. We also compare our predictions of 7-ray fluxes with the observed 7-ray 
light curves. In section 4 we discuss the implications of this work in constraining and 
understanding the dredge-up. 

2 Determination of 56 Ni dredge-up in SN 1987A us- 
ing the He I 10830 A line 

In spectra of SN 1987A obtained in the period 76-135 days post-explosion, Elias et al. 
(1988) and Meikle et al. (1989) identified a strong absorption feature at ~10,700 A with 
the blue-shifted P Cygni trough of He I 10830 A. Using the spectral synthesis model 
described in Fassia et al. (1998), where the lines are formed by pure scattering above 
the photosphere and where line blending is taken into account, we have confirmed this 
identification and shown that the He I line was blended with C I 10695 A. The He I line 
trough was blueshifted by about 5000 km/s. This is in sharp contrast to other infrared 
P Cygni lines where velocities of less than 2000 km/s were observed (Elias et al. 1988, 
Meikle et al. 1989). The high blueshift and depth of the trough indicate that the line 
was optically thick, and formed in the outermost layers of the expanding supernova 
envelope. 

As explained above, the presence of the He I line implies re-ionisation by the 7-rays 
of 56 Co. Modelling of the He I line profile and evolution was carried out as follows. 
We began with the unmixed explosion model 10H (Woosley 1988), calculated the ra- 
dioactive energy deposited in the outer envelope of the supernova and hence found 
the ionisation balance as a function of radius. We then determined the 2s 3 S popula- 
tion density profile. This was fed into our spectral synthesis model and the synthetic 
line profile was compared with the observations. As will be shown below, model 10H 
drastically underproduced the He I 10830 A feature, suggesting that upward mixing of 
radioactive material must have occurred during the explosion. To determine the degree 
of this dredge-up we gradually increased the mixing of the 56 Co until a match with the 
observations was obtained. 



2.1 Calculation of the energy deposition and the emergent 
luminosity of the 7-rays 

7-ray energy deposition was computed by performing a radiation transport solution 
for each 7-ray line produced in the decay 56 Co— > 56 Fe. Line energies and emission 
probabilities were obtained from Lederer & Shirley (1978). We also included 7-rays 
from the annihilation of positrons emitted in 19% of the 56 Co decays. For this reason, 
following Leising & Share (1990), we used a branching ratio of 0.38 for 7-rays with 
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energies of 0.511 MeV. 



We estimated the 7-ray energy deposition and the escape fractions, using a purely 
absorptive transfer equation. While Monte Carlo simulations provide the most realistic 
and physically accurate description of the 7-ray energy deposition, they are excessively 
demanding of computational resources. In contrast, the radiative transfer method, 
adopted in this study, is computationally efficient and has been shown to give results 
reasonably close to those obtained by the Monte Carlo techniques (Swartz et al. 1995). 
This is because the 7-ray energies are such that Compton scattering dominates. The 
ratio of scattered to incident photon energies is given by 



E v 1 + (£„/m e c 2 )(l -cos0) v ' 

where E v is the initial photon energy, and E' v is photon energy after it scatters from an 
electron at rest through angle 0. Thus the energy transfered in a single scatter can be 
large e.g. a ~2 MeV 7-ray photon gives up >0.8 of its energy when it scatters through 
> 90° (Swartz et al. 1995). Even though the scattering cross section is forward 
peaked, forward-scattering transfers very little energy to the gas. Therefore, one can 
think of the 7-ray line photon as travelling along a linear path until it suffers a large- 
angle scatter, at which point it gives up most of its energy and becomes a "continuum" 
photon. Consequently, in calculating the line transfer, the Compton scattering can be 
treated as producing a purely absorptive opacity. Thus, the radiative transfer equation 
is linear for the 7-ray lines. We can therefore calculate the energy deposition and 
escape intensities by direct integration of the radiative transfer equation. The emissivity 
for this calculation is given by the local volume rate of 7-ray production. Following 
Woosley's et al. (1994) approach, the opacity is assumed to consist of scattering and 
absorptive parts. The scattering component is 



N e a KN (E u ) 



scat l + (^/m e c 2 )(l-cos^ cat ) 
and the absorptive part is: 



(2) 



where E v is the initial photon energy, N e is the total (bound and free) electron density, 
<jkn{E u ) is the Klein-Nishina cross section, and scat is the scattering angle. In our 
calculations we have followed (Woosley et al. 1994) in assuming that all photons scatter 
through the same angle, 9 scat =90 o . This assumption gives results that agree to within 
a few percent with those obtained from more detailed Monte-Carlo calculations. 

The code solves the frequency dependent equations of radiative transfer by integrating 
along impact parameters parallel to the observer's line of sight (see Figure 1). Following 
Swartz et al. (1995) the transfer equation in the "(p,z) representation" along a ray is: 

f)T ± 

±-£-(p,z)=r h ,- Kv (r)I? (4) 
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where I + (I~) denotes the outgoing (ingoing) radiation, z is the position along the ray 
(Figure 1) and r\ v = f rad Sj(E v )p/4ir is the local 7-ray emissivity. Here, f ra d is the initial 
mass fraction of 56 Ni, p is the mass density and s^(E v ) is the time-dependent rate of 
energy release per gram of radioactive material: 

1.25xl0 17 E v b Ev e^t , _! _ x 
^(Eu) = (erg g rad s ) (5) 

r C'o 

where b Ev is the branching ratio for a line of energy E v , and rc is the mean life time 
of 56 Co (rc =113.7 days). Introducing the optical depth along the ray, dr = —K u dz, 
the equation of transfer becomes: 

dr k v 

If we now take into account that iV e = (p/m u )Y e where m u is the atomic mass unit and 
Y e is the total number of electrons per baryon, we have: 

d1 v _ frad ~-± ,„\ 

dr ~ Y e +± [n 

where 

I=(^-)L (8) 

s 7 m u 7V e 

Since no radiation is incident from outside the ejecta the boundary conditions are: 7 _ =0 
at the upper boundary and from symmetry at z — 0, the lower boundary condition is 
I~(p,z) = I + (p,z). Integration of equation (8) gives: 

Jfa) = X(r m )e- A - + T +1 l^e-^dt (9) 

where Tj and r i+1 are the optical depths to points i, i+1 along the ray and Ar=(r i+1 — Tj). 
The impact parameter grid consists of rays tangent to concentric shells (Figure 1), with 
a single ray passing through the center. Since f ra d, Y e are constant within each shell, 
equation (9) can be integrated exactly. 

By setting k u = K v abs we can calculate the rate (Re) at which the energy is deposited 
locally (in erg cm -3 s" 1 ) from: 

Re = ^ttk v J v (10) 

v 

where J u is the mean intensity of the 7-rays with initial energy E v . The emergent 
luminosity of the individual (unscattered) 7-ray line photons can also be calculated 
with the 7-ray transport calculation described above by setting k v = N e x a KN (E v ) 
(Woosley et al. 1994). 
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2.2 Calculation of the He I 10830 A line 



In a H/He envelope of low ionisation, all the energy released by the Compton scattering 
of the 7-rays is channelled into ionisation of the H and He, rather than into heating 
the electron gas (Meyerott 1980). Thus, we assume the key populating process for the 
helium levels is recombination. We ignore direct excitation by fast electrons since they 
only excite singlet states, which decay rapidly back to the ground state. Of the recom- 
binations to excited levels of He I, approximately three-quarters are to triplet states, 
with the remainder going to singlet states (Osterbrock 1989). Recombinations to singlet 
states cascade rapidly to the ground state. However, because there is a substantial op- 
tical depth in the ls 1 S-np 1 P transitions (Graham 1988) all singlet recombinations will 
eventually pass through the 2s 1 S state. Atoms in the 2s 1 S level decay by two-photon 
emission (A=51 s _1 ) to the ground state. Therefore, only a small population of excited 
singlet states is built up. 

Recombinations to triplet states lead, through downward radiative transitions, to the 
highly metastable 2s 3 S level where the population can become quite substantial. As 
described in Fassia et al. (1998), there are a number of mechanisms which can depop- 
ulate this level. Firstly, a very weak single-photon radiative decay to ls 1 S can occur 
(A = 1.27 x 1CT 4 s^ 1 ) (Osterbrock 1989). Of considerably greater significance, however, 
is the 2-stage intersystem radiative decay 2s 3 S^2p 3 P^ls 1 S. This occurs in the pres- 
ence of a radiation field (such as from the photosphere), which is required to excite the 
first stage. There are also two important depopulating processes involving collisions. 
One of these involves thermal electron collisions causing excitation or de-excitation 
from 2s 3 S across to singlet states. The other process is Penning ionisation (Bell 1970, 
Chugai 1991): 

He I (2 3 S) + H -»■ He I (l 1 ^) + H + + e. (11) 

which takes place when hydrogen and helium are microscopically mixed. Therefore, 
the population balance of the 2s 3 S level is described by the equation: 

n He (2»s)(n e Q + C P + R + A) = a(n 3 L)n e n He n- (12) 

where Q is the sum of the collision rates from 2s 3 S to all singlet states, and has the value 
1.826 xl0~ 8 cm 3 s _1 (Berrington & Kingston 1987). Cp is the Penning ionisation rate 
and is given by C P = ^pUu where 7p = 7.5 x lO -1 ^^^) 1 / 2 cm 3 s _1 (Bell 1970). For the 
recombination temperature of hydrogen (T~5000 K), which we adopted for the days 
concerned (Catchpole et al. (1987)), 7p=3x 10~ 9 cm 3 s _1 . a(n 3 L) = 3.26x 10 -13 cm 3 s -1 
is the total recombination coefficient for the triplet states at a temperature of 5000 K 
(Osterbrock 1989). A is the spontaneous transition probability and R is the two stage 
inter-system radiative decay rate (Chugai 1991). This is given by: 

R(v) = B 23 ^ t yF?e°^A^A 31 p 31 (13) 

where (3 3 i is the escape probability for the photon emitted in the decay 2p 3 P^ ls 1 S and 
is the flux of the continuum at A=10830 A, determined from the observed infrared 
spectra. A(A) is the extinction calculated using the empirical formula of Cardelli et al. 
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(1989) assuming Ay =0.6 (Blanco et al. 1987). For the distance D, a value of 50 kpc 
was adopted. 



For a deposition rate e» (erg cm 3 s r ) in species i, the abundance of the next ionisation 
stage i+1 is given by the energy balance : 

n i+1 = — ^— (14) 
n e aiWi 

where n e is the electron density, Oii is the recombination rate to all levels of the species 
i and Wi is the energy required to produced an ion-electron pair. Considering an enve- 
lope that consists only of hydrogen and helium and taking the Penning ionisation into 
account, equation (14) reduces to: 

1- n He t 2 3 S) n H jC P = n HII n e a H j (15) 

w H 

= n He n n e a He (16) 

where e = e# e + eh is the energy deposited in the envelope. The ratio of the energy 
deposited in H and He is e He /^H = rnY , where Y = n He /n H . From the Bethe-Bloch 
formula for energy loss of fast electrons in H/He material, we find m to be 1.7. Assuming 
n e = n>Heii + n<Hii we can then calculate the population of the 2s 3 S level, n 2s 3 5 (R) as 
a function of radius, R, using equations (15), (16) and (12). This requires the energy 
deposition rate e(R) which we obtain from the radiative transport calculation described 
in 2.1, the total number density n to t(R) = n He i(R) +nn i(R) +n e (R), and the relative 
abundance of hydrogen and helium, Y = n Heii(R)+n H ei(R) _ rp^g \ &s ^ ^ wo q Uan tities were 

J to ' n HII {R)+n HI {R) ^ 

obtained from detailed explosion models. 

The calculated density profile of the population of the 2s 3 S level, n 2s 3 5 (R), was then 
fed into our spectral synthesis code (Fassia et al. 1998), thus producing a model He I 
10830 A P Cygni line profile. This was then compared with the observed profile. 



3 Results 



We first applied the technique described above to the radially-symmetric, unmixed 
model 10H (Woosley 1988). This model represents the explosion of a blue supergiant 
consisting of a 6 M Q helium core, and a 10 M hydrogen envelope. The model incor- 
porated an explosion energy of 1.4 x 10 51 ergs and a 56 Ni mass of 0.07 M . Among 
the unmixed explosion models, 10H was one of the most successful at reproducing the 
observed UVOIR light curve. Nevertheless, significant discrepancies existed. Model 
10H also failed to explain the early detection of X-rays and 7-rays. Nor did it account 
for the high velocities observed in the iron-group line profiles. The present work shows 
(Figure 2a) that model 10H also fails to reproduce the pronounced He I absorption 
trough observed at ~10700 A. The lack of this He I trough means that the emission 
component of the C I 10695 A line is unsuppressed resulting in a poor match to the 
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observed spectrum. This discrepancy is not surprising since, in 10H, Ni only extends 
to 1100 km/s (Figure 3), whereas at 76 days for example, the helium photosphere is at 
~ 4700 km/s. The 7-rays simply do not penetrate to sufficiently high velocities. 

In order to account for the early detection of the X-rays and 7-rays Pinto & Woosley 
(1988) introduced ad hoc outward microscopic mixing into model 10H, yielding model 
10HMM. In this model, radioactive material is mixed outwards through the helium core 
and into the hydrogen envelope. Model 10HMM successfully reproduced the UVOIR 
light curve and accounted for the early and prolonged detection of X-rays and 7-rays. 
However, on applying our technique to this model, we found that it provides a poor 
match to the observed He I 10830 A line profile (Figure 2b). This is because the 
extensive dredge-up of 56 Ni invoked in model 10HMM (Figure 3) produces a 2s 3 S 
population density that is much too high in the line forming region. 

In order to attempt to match the observed line profile we repeated the microscopic 
mixing procedure that Pinto & Woosley applied to model 10H (Pinto, private com- 
munication). Turbulent or macroscopic mixing will be discussed in the next section. 
Starting from the centre, we specified a velocity interval and homogenized the material 
composition over that interval. We repeated this procedure, moving outwards until 
the edge of the ejecta was reached. This brings a fraction of the radioactive material 
to high velocities. Nevertheless, the abundance gradient is such that the bulk of the 
radioactive material remains at low velocities. 

For each of the mixed models produced we calculated the 7-ray energy deposition, ioni- 
sation balance and 2s 3 5 population density as a function of radius, and then compared 
the resulting synthetic He I 10830 A line profile with the observations. Figure 3 shows 
the 56 Ni distribution for two of our mixed models (denoted 10HMA and 10HMB) that 
provided a good match to the observed He I 10830 A line profile. The synthetic he- 
lium line profiles deduced using model 10HMB are compared with the observations in 
Figure 4. Equally good fits were obtained for model 10HMA. 

We also modelled the 7-ray line light curves for models 10HMA and 10HMB using the 
radiative transport calculation explained in 2.1 and adopting k v = N e x o"xat(-Eo)- 111 
Figure 5, these light curves are compared with 7-ray observations (Leising & Share 
1990; Sandie et al. 1988), assuming a distance of 50 kpc. Also shown are the 7-ray 
light curves deduced using model 10HMM. 



4 Discussion 

The results show that we can reproduce the He I line by reducing the amount of 
microscopic mixing that was invoked in model 10HMM. Indeed, in model 10HMA 
there was negligible radioactive material at velocities exceeding ~3,900 km/s. This is 
in contrast with 10HMM where ~1% of the total 56 Ni mass has velocities in excess 
of 3,900 km/s. We find better agreement with the more recent work of Nomoto et al. 
(1998). Based on their explosion model 14E13, their bolometric light curve calculations 
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suggest that 56 Ni is mixed up to 3,000-4,000 km/s. 

For models 10HMA and 10HMB we find that, respectively, 4% and 3% of the total 
56 Ni mass (0.07 M Q ) lay above 3,000 km/s. This is consistent with late-time infrared 
spectroscopic observations which showed that a small fraction of the cobalt and iron 
was travelling at velocities ~3,000 km/s (Spyromilio, Meikle & Allen 1990; Haas et al. 
1990; Tueller et al 1990). From their 18 fxm and 26 //m [Fe II] line profiles, Haas et al. 
inferred that > 4% of the iron mass had an expansion velocity greater than 3,000 km/s. 
We also note that 7-ray light curve predictions based on models 10HMA and 10HMB 
are consistent with the observations (Figure 5). 

Models 10HMM, 10HMA and 10HMB are equally successful at accounting for the 7- 
ray light curves and late-time infrared line velocities. However, in reproducing the 
observed early-time He I 10830 A line, models 10HMA and 10HMB are clearly superior 
to 10HMM. We conclude that, given current technology, the He I line technique provides 
a better measure of the dredge-up of radioactive materials in type II supernovae. Some 
of the advantage of the He I line method stems from the fact that it uses near-IR 
observations taken at relatively early times when the supernova is bright. An additional 
advantage is that, since the method is based primarily on the optical depth of the He I 
P Cygni line, it is insensitive to fluxing errors (Fassia et al. 1998). 

In the above work, microscopic mixing was assumed, and so Penning ionisation domi- 
nated the depopulation of the 2s 3 S level (c/ Fassia et al. 1998). However, hydrodynamic 
calculations show that mixing is turbulent and takes place on macroscopic length scales. 
Multi-dimensional studies predict that hydrogen and helium bubbles will be dragged 
towards the inner parts of the ejecta, while clumps of helium and heavier elements, 
including radioactive nickel, will penetrate the hydrogen envelope (Fryxell, Miiller and 
Arnett 1991, Herant & Benz 1991). Pure helium clumps in the hydrogen envelope will 
tend to dominate the population of the 2s 3 S level as they are not subject to the effects 
of Penning ionisation. Helium that is microscopically mixed with hydrogen will make 
a negligible contribution even if its abundance is much greater. In Fassia et al. (1998) 
we demonstrated that the presence of helium clumps in the hydrogen envelope reduces 
significantly the degree of 56 Ni dredge-up that is needed to drive the He I 10830 A line. 
In that study we assumed that the 56 Ni was microscopically mixed with the helium 
bubbles and the hydrogen envelope. However, in the case of SN 1987A, there is good 
evidence that some of the 56 Ni was ejected in fast-moving clumps (e.g. Spyromilio, 
Meikle & Allen 1990). Thus, if the 56 Ni in the vicinity of the He I line forming region is 
not microscopically mixed, it may be necessary to increase the degree of 56 Ni dredge- 
up to provide the required rate of helium ionisation. Clearly, this is a complicated 
problem, requiring the use of 2-D and 3-D explosion models. We therefore defer to a 
later paper an examination of the implications of macroscopic mixing in SN 1987A for 
the He I technique. This should ultimately allow us to place constraints on the initial 
inhomogeneities. 

Finally, we note that the He I observations take place well before the 7-ray lines reach 
maximum brightness. Thus, by observing the He I 10830 A line in type II supernovae 
at early times, we should be able to make accurate predictions of the 7-ray fluxes, 
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providing a valuable "early-warning" facility for 7-ray observers. 
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Figure Captions 



Figure 1 : The coordinate system (p,z) used for the radiative transfer calculation. The 
dashed line shows an example of an impact parameter ray along which the frequency 
dependent equations of transfer were integrated. 

Figure 2 : Comparison of observed spectra of SN 1987A at 105 days, with synthetic 
spectra whose He I 10830 A line profiles are based on models 10H and 10HMM. Data 
are from Elias et al. (1988). 

Figure 3 : Distributions of 56 Ni for models 10H (dashed line), 10HMM (dotted line), 
10HMA (solid line) and 10HMB (dot-dashed line). 

Figure 4 : Comparison of observed spectra with synthetic spectra whose He I 10830 A 
line profiles are based on model 10HMB. Data are from Elias et al. (1988) and Meikle 
et al. (1989). 

Figure 5 : Comparison of the 7-ray line light curves of model 10HMA (solid line) and 
model 10HMB (dotted line) with the observations. Data are from Leising & Share 
(1990) (open triangles) and from Sandie et al. (1988) (open squares). The dashed line 
represents the predictions of model 10HMM. 
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